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C , One-dimensional non-equilibrium kinetic Ising models evolving under the competing effect of spin 

• flips at zero temperature and nearest neighbour spin exchanges exhibiting a parity-conserving (PC) 

' ' ' phase transition on the level of kinks are investigated here numerically from the point of view of the 

0^ , underlying spin system. The dynamical persistency exponent and the exponent A characterising 

' the two-time autocorrelation function of the total magnetization under nonequilibrium conditions are 

(— I ' reported. It is found that the PC transition has strong effect: the process becomes non-Markovian 

. and the above exponents exhibit drastic changes as compared to the Glauber-Ising case. 
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, In recent years , two nonequilibrium dynamical critical exponents have been discovered, which arise under non- 
^2 ' equilibrium conditions. The non-equilibrium (short-time) exponent A characterizes two-time correlations in systems 
^ [ relaxing to their critical state in the process of quenching from infinitely high temperatures to Tc [^jj^] . Recently, one 
more critical exponent was proposed the persistence exponent 0, associated with the probability p(t) cx t~^, that 
the global order parameter has not changed sign up to time t after a quench to the critical point Q . For some known 
examples, cited in ref. H, the scaling law 
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eZ = A-d+l-| (I) 



^ • is satisfied (here d is the dimensionality and 77 is the static critical exponent of the order parameter correlation 
CO [ function), which has been derived assuming that the dynamics of the order parameter is a Markovian process. In 
25 ■ general, however, A and Q have been proposed Q| to be independent, new critical dynamical exponents. 

One of the soluble examples is the d — 1 Ising model with Glauber kinetics. In this case the critical temperature 
is at r = 0, and as shown in ref. the persistence exponent is = 1/4 for the global order parameter which is 
the total magnetisation M{t). Moreover, A is known to be A = f in this model. The aim of the present note is 
^\ to study these new dynamical critical exponents in a simple non- equilibrium Ising system (NEKIM) introduced in 
"Y^ I^.The phase diagram of NEKIM consists -essentially - of a line of (first order) Ising-type transitions, which line ends 
I at a specific point. At this endpoint, on the level of kinks (phase boundaries) a second order phase transition takes 
^ I place from an absorbing to an active state which, however, belongs to the parity conserving (PC) universality class 
I ' P,[7p|J^-p^. The critical fiuctuations of this PC transition exert a pronounced effect on the underlying spin system 
^5 , as found earlier |jll| thus e.g. the the classical dynamical exponent Z, defined, as usual through t oc with ^ cx p~'' 
Q ■ was found to be Z = 1.75(f) instead of the Glauber-Ising value of Z = 2. (We note here that p — e~'>^ plays the role 
Q of '^zF" in Id with Tc = and the static exponents are defined as powers of p for T ^ ). 



The question arises how the citical fluctuations of the PC transition affect the other two critical dynamical exponents 
Q and A. Before entering into the details of our results for @ and A, the model will be described in some detail. 



^ . In NEKIM the system evolves under a combined effect of spin-flips and spin-exchanges. The spin-flip transition 
rate in one-dimension for spin Si (s^ = ±1) sitting at site i is |f 2[| : 
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(f + Ssi^iSi+i) (l - ^Si(si_i -I- Si+i)j (2) 



where 7 — tanh2J/fcr (J denoting the coupling constant in the Ising Hamiltonian), F and d are further parameters. 

^(1 - S) and ^( 



At T = 7 = f and there are two independent non-zero rates ^(1 — S) and ^(f + S), responsible for random walk. 



and pairwise annihilation of kinks, respectively. 
The spin-exchange transition rate of nearest neighbour spins (the Kawasaki |f 3| rate at T = 00) is 
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Wii+i = ^Pex[i — SiSi+i], where pex is the probabiHty of spin exchange. Spin-flip and spin-exchange have been apphed 
alternatingly. 

In this system, at T = 0, a PC type phase transition takes place. In j5| we have started from a random initial state 
and determined the phase boundary in the {S,pex) plane. In the following we will choose a typical point on this phase 
diagram and make simulations at this point. The parameters chosen are: F — .35,pex — .3, (5c = —.395(2). In the 
simulations the spin-flip part has been applied using two-sublattice updating. After that we have stored the states of 
the spins and made L (L is the size of the system) random attempts of exchange using always the stored situation for 
the states of the spins before updating. All these together has been counted as one time-step of updating. ( Usual 
MC update in this last step enhances the effect of Pex and leads to 6c — —.362(1)). 

In M it has been argued that for studying non-equilibrium critical dynamics, the global, rather than the local order 
parameter should be considered. The non-equilibrium nature of the problem under consideration is partly due to the 
model itself and partly due to the conditions of a quench from T = cxd to T = 0. This latter means that we will 
restrict ourselves to completely random initial states and follow the behaviour of the system using the rules described 
above. The persistency exponent is defined via the probability p{t) that the global order parameter, which in our 
case is the total magnetization: < Mk=o{t) >= < J2i ^ii^) >i ^i^s not changed sign up to time t: 

p{t) oc r'^ (3) 

Finite size scaling (FSS) applied to the persistence problem Q leads to the form: 

p{t) = L-'^'^git/L^) (4) 
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FIG. 1. p{t)L'', a = QZ, plotted against t/L^ with Z = 1.75 and 9 = .65, .67, .69. For clarity the .65 and .69 data 
have been multiplied and divided by a factor of 2, respectively. The simulational data exhibited here have been obtained for 
L = 100, 200, 400, 800 sized systems labeled by different symbols, for statistical averages between 3 x 10® — 10* samples. 



Simulations have been carried out in the range 50 < L < 2000 with periodic as well as antiperiodic boundary 
conditions and at least for 10^ independent runs. Fig.l shows those of our results which have the best statistics with 
averages over up to 10^ independent random initial configurations ( those configurations, however, for which M(0) = 
exactly, were discarded) with periodic boundary conditions for L — 100,200,400 and 800. Using Z — 1.75(1) jlj] the 
best fit corresponds to 8 = .67(1). The fact that periodic boundary conditions allow only an even number of kinks 
leads eventually to perfect ordering of spins seen as levelling off of curves on Fig. 1 . Our simulations with antiperiodic 
boundary conditions have led to the same value of Q as above though the form of the scaling function g{t/L^) in this 
case is different, of course. For comparison we have also simulated the exactly soluble Glauber-Ising case and found 
the expected value of = .25 within the accuracy of the simulations. 



The local autocorrelation function defines the new exponent A 




A{t,0) = -<J2s^{0)s^it)>xt-i. (5) 
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We have made simulations for this quantity after starting with a random initial configuration, and allowing the 
system to evolve according to the rule of NEKIM as described above. Averaging has been taken over random initial 
configurations in a chain of length L = 1000. The result is shown on Fig. 2. The best fit has been obtained with 
A = 1.49(3), using Z = 1.75(1). For comparison, numerical results for the corresponding quantity in the Glauber-Ising 
limiting case are also displayed on Fig. 2. It is worth mentioning, that data for t < 10 had to be discarded in both 
cases; power law behaviour is seen only for later times and this fact does not change if the number of averages taken 
is increased even by an order of magnitude. 




FIG. 2. Time dependence of the local autocorrelation function A{t, 0). L = 1000 and the number of averages over independent 
random initial states was 1.5 x 10^ also for the Glauber-Ising case which is shown for comparison. 

Following Q , we will study now the two-time autocorrelation function for the global order parameter: 
A^^°^°-'- {ti,t2) = L < Mk=o{ti)Mk=oit2) > or rather its normalized form, namely 

a{h,t2) - Aa'''''-\h,t2)/^/S(0^VS(^t2) = fir) (6) 

Here S{0,t) = L < -^Ei ■Si(i)]^ > is the structure factor at the ferromagnetic peak and the second equality follows 
from scaling assumption (see later in more detail). Morerover, for y — > oo, f{y) ~ y-i^-d+^-^/^)/^ is the expected 
power law behaviour. Nevertheless, if the process is Markovian, the power law behaviour of /(|^) has to hold for all 
ti > t2 as shown in [Q. 

The second moment of the global magnetization (structure factor) should behave as |l5) 

5(0,i)^t('^-^)/^^i^ (7) 

We have found earlier, in |jri| that /3 = .00(2), i.e. even at the PC point the Ising phase transition is of first order and 
thus 5(0, t) oc t^^^ . Moreover, via the above applied scaling law d — 2 + rj = 2^, rj = 1.0(1) follows at the PC point, 
too. 

Fig. 3. shows a{ti/t2) as a function of ti/t2 for six different values of ^2 = 3, 5, 10, 32, 50, 100. We have simulated 
yl^'°''"'(ti, t2) while for the denominator we have used the power law behaviour as indicated above with 1/Z = .57. 
Unfortunately, it is very hard to get Af'°'"^'(ii, t2) to a satisfactory accuracy. 
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FIG. 3. The normalized autocorrelation function a{t\/t2) at the PC point for six different values of t2 increasing in the 
downward direction in the range 3-100. The Galuber case is exhibited again for comparison. L — 1000 and the number of 
averages over independent initial states: 1 — 2 x 10"" . 

d-X 

The reason is probably the fact that in the scahng form of the leading order term is oc (|^) ^ which 

is non-singular in the present model ( it is marginal with d — X for the Id Glauber model while it is singular for the 
2d Ising case). This means that the k — mode is not special for the two-time structure factor (while it is for the 
equal-time structure factor) and thus a power law behaviour is correction to scaling p8[ |. In more detail: the scaling 
form for the two-time structure factor can be written as 

< M_fe(ti)Affe(i2) >= 

where L{t) ^ t^/^ and for fc — > fi{kL{t)) const. Moreover since the second term is correction to scaling, 
Ai < d ~ \. For ti — t2 we get the usual structure factor and for fc ^ without the correction term this is the 
form cited in below their eq.(16). In the present case, knowing that the singular term is missing we can say 
limfc^o fi{kL{t)) = 0, scaling can still be present and Ai plays the role of d — A (and we will use d — A for Ai in the 
following even for X > d). 

Now turning back to our Fig. 3, apart from the first three decades in time, fluctuations hinder drawing any con- 
sequence concerning the (correction to scaling ) behaviour of A^''°^°'''{ti,t2) even for averages of order 10^. For the 
quantity ti/t2 this fact narrows down the interval of analyzable data even more. Nevertheless, it is clearly seen 
that the dynamic scaling assumption expressed in eq.(6) ( and which eventually can be expected to hold only for 
ti ^ l,t2 ^ l!) starts to be fulfilled to an accuracy below 1% only for values ti > ^2 ^ 50. Actually this is not 
typical: it has been proposed [ p^ that in case of systems quenched to their critical temperature (and here Tc = 0) 
universality and scaling may appear in a quite early stage of time evolution, far from equilibrium, where ^(i) '~ t^^^ 
is still small. Based on the scahng relation for such early time intervals, a new way for measuring static and dynamic 
exponents has been proposed [p"5[p^ and applied also for the local autocorrelation function |0 . Some of our earlier 
results also show that, indeed, power law behaviour sets in for quite early times already. Thus e.g. in [ pl| concerning 
the structure factor S{0,t) ^ , the power law behaviour was apparent already for very early times and for such 
low values oi L as L — 128, provided the number of averages in the simulation was high enough (above 10^). The 
obtained result, x = .57 = 1/Z was made use of above. 
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h{kL{h))+t2' 



f2ikL{ti)), ii,i2>l, d>X 



(8) 
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For the sake of comparison we have carried out similar simulations of ^^'"''"'(^i, ^2) for the exactly soluble Ising- 
Glauber case {pex = 0: <5 = 0), some of these are also exhibited on Fig. 3. Here dynamic scaling is fulfilled (to similar 
accuracy as above) already for ti > t2 ^ 5 and the expected power law behaviour is seen within error. It is worth 
mentioning that similar value for A, i.e. A = 1.0 results from simulations in the whole absorbing region (thus e.g. for 
Pex = -35,(5 = 0). 
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FIG. 4. The dependence of the normalized autocorrelation function on ti/t2 in its first decade for t2 = 100. L = 1000 and 
the number of averages was lO'^. The straight line is best power law fit with (A — ^)/Z = .581 

In order to establish whether the process is Markovian or not at the PC point it will be sufficient to examine the 
first decade in the variable |^ in a region where dynamic scaling holds. Fig. 4 shows the result for the case t2 = 100 

and for averages over 10^ independent initial states, again taking L = 1000. For the exponent '■^'^^^^^ the value .58 
results as best fit, which, according to eq.(l) should equal Q. From here one arrives at A = 1.51(1) which is in accord, 
within error, with the value obtained above from the local autocorrelation function. Thus, supposing the Markovian 
property to hold has led to contradiction because the measured value of 8 is .67(1). 

The results together with critical exponents obtained earlier in pT| are summarized in Table I. 





/3 


7 




Z 


e 


A 


Glauber-Ising 





1/2 


1/2 


2 


1/4 


1 


PC 
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TABLE I. Simulation data for static and dynamic critical exponents for NEKIM 
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In summary, wc have carried out mimerical simulations to investigate the non-equihbrium dynamic critical exponents 
O and A with the aim to check the Markovian nature of the nonequilibrium Ising system in Id at the parity conserving 
phase transition point of the phase diagram of NEKIM. On the basis of of the present results we have been led to 
the conclusion that the effect of fluctuations felt by the spin system at the PC transition is such that the dynamical 
process becomes non- Markovian. The difference is quite pronounced, definitely beyond numerical errors. 
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